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ABSTRACT 

We report evidence of ordered orbital motion in luminous star-forming galaxies at 
z ^ 1.3. We present integral field spectroscopy (IFS) observations, performed with the 
OH Suppressing InfraRed Imaging Spectrograph (OSIRIS) system, assisted by laser 
guide star adaptive optics on the Keck telescope, of 13 star-forming galaxies selected 
from the WiggleZ Dark Energy Survey. Selected via ultraviolet and [Oil] emission, 
the large volume of the WiggleZ survey allows the selection of sources which have 
comparable intrinsic luminosity and stellar mass to IFS samples at z > 2. Multiple 
1-2 kpc size sub-components of emission, or 'clumps', are detected within the Ha 
spatial emission which extends over 6-10 kpc in 4 galaxies, resolved compact emission 
(r < 3 kpc) is detected in 5 galaxies, and extended regions of Ha emission are 
observed in the remaining 4 galaxies. We discuss these data in the context of different 
snapshots in a merger sequence and/or the evolutionary stages of coalescence of 
star-forming regions in an unstable disk. We find evidence of ordered orbital motion 
in galaxies as expected from disk models and the highest values of velocity dispersion 
(tr > 100 km s^^) in the most compact sources. This unique data set reveals that the 
most luminous star-forming galaxies at z > 1 are gaseous unstable disks indicating 
that a different mode of star formation could be feeding gas to galaxies at z > 1 , and 
lending support to theories of cold dense gas flows from the intergalactic medium. 

Key words: galaxies: high-redshift - galaxies: kinematics and dynamics - galaxies: 
formation - galaxies: evolution. 
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1 INTRODUCTION 

An open question remains regarding the accretion of gas and 
build up of stars in massive galaxies at the peak epoch of 
cosmic star formation. The traditional theories of isolated 
disk-formation res ulting from gas accretion from the halo 
IjEggen et al.lll962l l have been challenged by recent observa- 
tions and theories at high redshift. In new models of disk 
formation at high redshift, c ooling flows supply gas directly 
to the centres of galaxies (jKeres et alj l2005l : iDekel et alj 
l2009al ). Direct observations of the gas motions in gala:xies 
are crucial to test these new models and ultimately to dis- 
entangle how galaxies form at high redshift. 

Morphological studies made possible by Hubble Ad- 
vanced Camera for Surveys (ACS), Wide Field and Plan- 
etary Camera (WFPC) and other deep imagers have pro- 
vided the first clues to the processes governing high-redshift 
galaxy formation. High-redshift galaxies exhibit different 
morphologies to local galaxies, wit h higher fractions of both 
compact and irregu lar gala:xies (jGlazebrook et al.l 1 19951 : 
[Abraham et al.lll996r ). This deviation from the Hubble Se- 
quence indicates that different mechanisms may be driving 
star formation and galaxy evolution at earlier cosmic times. 
Initial interpretations of the complex structures seen in the 
Hubble Ultra Deep Field (HUDF) favor higher rates of ma- 
jor mergers, a natural consequence of hierarchical merging 
models in a ACDM Universe ( e .g., i Consclicc ct al. 2003; 
iLotz et al.l l2006l : lOverzier et all [2OO8I : lEngel et all I2OI0I ). 
Yet, despite large dense regions in these galaxies, observa- 
tions of their stellar populations r eveal galactic properties 
reminiscent of local disk galaxies ( Elmeg reen fc Elmegreenl 
I2OO5). In this case the large dense regions or 'clumps' of 
star formation may be related to HII regions found in 
spiral galaxies but o f much larger size (jJones et al.l l20ld : 
ISwinbank et al.l [20101 ). Models have been developed to ex- 
plain the clumpy structures in disks and their evolution 
into galactic structures seen at 2 = (e. g., Noguchil 19991: 



Immeli et ai]|2004al : iBournaud et al.ll2007l : lElmegreen et al.1 
20081 ). In these models, turbulent gas collapses from Jeans 



instabilities resulting in massive star-clusters that, on a 
timescale of 1 Gyr, coalesce to the centre of the galaxy to 
form a pseudo-bulge. 

Morphological studies alone do not provide enough in- 
formation to distinguish between the two kinematically dis- 
tinct formation theories of merging and isolated disk forma- 
tion at these redshifts. Given the complex structure of the 
galaxies it is difficult to derive the major axis for slit po- 
sitioning making spectroscopic obse rvations less effective i n 
discerning rotation signatures fe.g. lErb et al.ll2004l . 120061 ). 
Spatially-resolved spectroscopy is therefore essential for dis- 
entangling kinematics and star-formation histories of galax- 
ies at high redshift. 

The past five years have seen a great number of inte- 
gral field spectroscopy (IFS) studies of the neb ular emission 
line, Hg, in g > 2 star-forming g alaxies ( e.g.. ICenzel et al.l 
20061. I2OO8I: iForster Schreiber et a l. 2009 : JLaw et al.ll2009l: 



Wright et al.ll2007l. 120091: lEpinat et al.ll2009l:lNesvadba et all 



20081 ) . These redshifts are favored due to the large number of 
Lyman break galaxies (LBGs) and sub-millimeter galaxies 
(SMGs) being discovered by methods tailored for high red- 
shift as well as the accessibility of Ha in the near-infrared 
(NIR). 



A main focus of these studies is the classification of 
galaxies kinematically to define clear observational distinc- 
tions between disks and interacting systems. Mergers ob- 
served with integra l field units fIFUs: iNesvadba et al. 20081 : 
IWright et al.ll2009l ) often show clear kinematic steps across 
one resolution element, indicating separation between two 
kinematically separate systems, whilst disks are identified 
by rotational signatures found in both velocity and velocity 
dispersion maps. A growing number of massive disks have 
been found with IFS at 2: > 2 which are well fit by disk 
models and which have Toomre parameters of Q > 1 indi- 
cating stable rotationally supported disks, lending observa- 
tional support to theories of rapid accre tion of cold gas and 
secular disk evolu t ion at high redshift (ICenzel et al.l 20081 : 
lEpinat et al.ll2009l : iLemoine-BusseroUe fc Lamareillell20ld ). 

However, classifying galaxies in these two distinct 
groups is an over-simplification. Many of the samples present 
a variety of kinematic types including dispersion-dominated 
systems, stable disks, unst able disks, minor mergers, merg- 
ers, and merger remnants. iLaw et al.l ((20091) present a sam- 
ple of lower mass (M*,at)g = 2 x 10^" M©) galaxies with kine- 
matics dominated by dispersion. Their results can be inter- 
preted as merger remnants, however they also caution that 
this is statistically unlikely and find that the kinematic clas- 
sifications are a crude over-simplification. They argue that 
their galaxies are observed at a range of times during their 
formation, and that their kinematic signatures are a result of 
how they acquire their gas. Despite the different kinematic 
interpretations, all high-redshift IFS samples show high dis- 
persions (a > 90 km s~^) with many (20—100%) systems 
showing no ordered rotation making classifications uncer- 
tain (e.g. I Cenzel et al.1 I2OO8I: iForster Schreiber et all l2009l : 



iLaw et al.ll2009l : IWright et al.ll2009D . 



With observation times ranging from 2000s - 40000s, it 
is only possible to follow up a small fraction of the bright- 
est galaxies at these redshifts. To supplement high-redshift 
kinematic results other studies target lensed galaxies mag- 
nified in both size and flu x, allowing the study of sub-L« 
galaxies (IJones et al.ll2010l) . Using this strategy, IFUs have 
probed all the way ou t to z = 4.9 with NIFS on Cemini 
(|Swinbank et al.l 120091 ) and have targeted stro ngly lensed 
galajc ies aX z = 1.4 — 3 . 1 with OSIRIS on Keck (jStark et al.l 
20l3g|: |jones et al.ll2010l : rYuan et al.ll201lh . With galaxy mag- 
nification individual clumps have been resolved with diame- 
ters of ~500 pc which will be essential in testing the models 
of clumpy d isk formation (jJones et al.ll2010l : ISwinbank et al.l 
|2009| . l20ld ). However, more kinematically resolved clumps 
are needed than are available from current studies (~ 20 at 
2: > 1) in order to rigorously compare with models. 

Despite great advantages, lensed studies require strong 
lenses with confirmed redshifts which can be observationally 
expensive in addition to requiring reliable lensing models 
and reconstruction. Furthermore, rather than bridge the gap 
between the z ~ 2 — 3 results and local galaxies, the majority 
of current lensing studies push out to higher redshifts. 

In order to correctly interpret IFS data at high red- 
shift it is critical to utilize low-redshift IFS results. By com- 
paring directly to these data and by artificially redshifting 
them out to z > 1 a better understanding of different kine- 
matic signatures in the high-redshift data can be reached. 
Recent efforts by lEpinat et al.l (|20ld ) have supplied a siz- 
able sample of FabryPerot 3D data of Ha emission in 153 
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Table 1. OSIRIS observations of WigglcZ galaxies 



WiggleZ ID 


ID 


R.A. 


Dec. 


z" 


Obs. 


^cxp. 


Exp. 


Filter 


Scale 






(J2000.0) 


(J2000.0) 




Date 


(s) 


Seq. 




(mas) 


Detections 


R03J032450240- 13550943 


WK0912_13R 


03:24:50.240 


-13:55:09.429 


1.288 


2009/12 


3600 


4 X 


900 


Hnl 


50 


S15J145355248-00320351 


WK1002_61S 


14:53:55.247 


-00:32:03.501 


1.305 


2010/02 


5400 


6 X 


900 


Hnl 


50 


R01J005822757-03034040 


WK0909.02R 


00:58:22.757 


-03:03:40.400 


1.362 


2009/09 


5400 


9 X 


600 


Hn2 


50 


R03J032206214- 15443471 


WK0909.06R 


03:22:06.214 


-15:44:34.711 


1.461 


2009/09 


8100 


9 X 


900 


Hn3 


50 


S15J144102444+05480354 


WK0905.22S 


14:41:02.443 


05:48:03.540 


1.282 


2009/05 


3600 


4 X 


900 


Hnl 


50 


R00J232217805-05473356 


WK0912.01R 


23:22:17.805 


-05:47:33.560 


1.297 


2009/12 


3600 


4 X 


900 


Hnl 


50 


S09J091517481+00033557 


WK1002.41S 


09:15:17.483 


00:03:35.570 


1.334 


2010/02 


7200 


8 X 


900 


Hnl 


50 


S00J233338383-01040629 


WK0809.02S 


23:33:38.386 


-01:04:06.290 


1.452 


2008/09 


5400 


6 X 


900 


Hn3 


50 


S11J101757445-00244002 


WK1002.14S 


10:17:57.444 


-00:24:40.020 


1.305 


2010/02 


5400 


6 X 


900 


Hnl 


50 


S15J142538641+00483135 


WK1002.18S 


14:25:38.639 


00:48:31.350 


1.305 


2010/02 


2700 


3 X 


900 


Hnl 


50 


S09J090933680+01074587 


WK1002.46S 


09:09:33.680 


01:07:45.870 


1.308 


2010/02 


3000 


5 X 


600 


Hnl 


50 


S11J110405504+01185565 


WK1002_52S 


11:04:05.504 


01:18:55.650 


1.320 


2010/02 


7200 


8 X 


900 


Hnl 


50 


S09J090312056-00273273 


WK0912_16S 


09:03:12.056 


-00:27:32.730 


1.327 


2009/12 


8100 


9 X 


900 


Hnl 


50 


R00J233618953- 10580749'' 


WK0809.05R 


23:36:18.955 


-10:58:07.321 


1.402 


2008/09 


2700 


3 X 


900 


Hn2 


100 


R00J232248993- 12151943^ 


WK0809.04R 


23:22:48.992 


-12:15:19.439 


1.455 


2008/09 


2700 


3 X 


900 


Hn3 


50 


Mis-Identifications 


R22J214700239+01333423 


WK0809.01R 


21:47:00.239 


01:33:34.230 


1.466 


09/2008 


3600 


4 X 


900 


Hn3 


50 


S00J234420801-09322580 


WK0809.03S 


23:44:20.801 


-09:32:25.800 


1.488 


09/2008 


2700 


3 X 


900 


Hn3 


50 



"WiggleZ Dark Energy Survey spectroscopic redshift. 

''Ho detected but observations unusable due to poor AO correction as a result of stray light on the low bandwidth wavefront sensor. 



nearby isolated galaxies. Disk models are fit to the local 
population and to the same galaxies artificially redshifted to 
z = 1.7, providing guidance in the best disk models to use 
to fi t to high-redshift kinem a tic re sults. I Green et al.l (|2010l ) 
and [Rodriguez- Zaurin et al.l (|201in also provide local com- 
parison samples for Ha velocity dispersion and morphology. 
Properties observed locally at parsec to kiloparsec resolu- 
tion are crucial for the interpretation of high-redshift IFU 
observations. 

Connecting the z < 0.2 and z > 1.5 results is a rela- 
tive dearth of IFS studies closer to 2 ~ 1. At this redshift, 
galaxies are at the tail end of the epoch of mass assembly 
but at more accessible distances for current instrumentation 
(where (1 + z)* surface brightness dimming is less severe). 
Yet suitable galaxies for such studies are rare as once the 
main epoch of massive galaxy formation is over at 2: ~ 2, 
the abundance of galaxies vv^ith high star formation rates 
falls rapidly. They are difficult to select with the Lyman- 
break technique at z ~ 1 as the drop in emission blueward 
of 912A cannot be detected with the low sensitivity and res- 
olution of ultraviolet (UV) instruments. These galaxies can, 
however, be found in very large high-redshift surveys, such 
as the WiggleZ Dark Energy Survey. The WiggleZ survey 
contains a rare population of UV-luminous galaxies out to 
z = 1.5 with star formation rates (from nebular emission) 
of ~100 Mq yr~^. These galaxies may be the link needed 
to bridge the gap between observations of local luminous 
galaxies, and LBGs and SMGs found at higher redshift. In 
this paper, we present integral-field spectroscopy of the Hq 
emission line in 13 star-forming galaxies from the WiggleZ 
Dark Energy Survey at z ~ 1.3. 

In this paper, we present 13 galaxies and their physical 
properties. In Section 2 we describe our sample selection, ob- 
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Figure 1. Nebular luminosity as a function of redshift for the 
largest high-redshift IFS surveys with the work presented here 
shown by the red diamonds. The gray points show the [OH] lu- 
minosity of other WiggleZ galaxies c onsidered for this s ample. 
The red circles show the sample from lEpinat et al.l 112003) taken 
from the VVDS survey and observed with SINFONI o r i the V LT. 
The green circles show the sample from iLaw et al.l 1120091) ob- 
served with OSIRIS on Keck. The blue inverted trian gles show 
the SINS sample from iForster Schreiber et al.l 1 120091) observed 
with SINFONI, and the green triangles show the sample from 
IWright etal] ||2009| ) observed with OSI RIS. The orange points 
are SMGs with measure d Ha taken from lSwinbank et al.l 1120041 ) 
and lTakata et al.l (l200d) . 
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servational technique, and data reduction methods. In Sec- 
tion 3 we discuss the properties of the galaxies, in particular 
we describe the morphological and kinematic properties. In 
Section 3.6 and 3.7 we discuss the methods of disk fitting 
and surface brightness profile fitting analysis. The individual 
and global properties of the 13 systems described in Section 
3 are analyzed in Section 4 in the context of theories of 
galaxy formation at high redshift and the implications for 
the mechanism triggering massive rates of star formation. 
We conclude with a summary of our results and interpre- 
tations in Section 5, giving a brief list of the key points 
presented in this paper. The properties of the 13 individual 
galaxies are described individually in the Appendix. 

A standard ACDM cosmology of fim = 0.3, r^A = 
0.7, h = 0.7 is adopted throughout this paper (e.g. 
iKomatsu et al.ll201lh . In this cosmology, at redshift z — 1.3, 
1 arcsec corresponds to 8.38 kpc. 



2 OBSERVATIONS 

2.1 Target Selection 

Our sample of 13 ga laxies was selected from the WiggleZ 
Dark Energy Survey (jDrinkwater et al.ll201Cl ). WiggleZ is a 
UV-selected survey of 238,770 unique emission line galax- 
ies in the redshift range < z < 1.5, measured over 245 
nights with the AAOmega spectrograph (jSharp et al.ll2006l ) 
on the 3.9-m Anglo- Australian Telescope (AAT). The pri- 
mary aim of the survey is to precisely measure the scale 
of baryon acoustic oscillations imprinted on the spatial dis- 
tribution of these galaxies at look-back times of 4-8 Gyrs 
IjBlake et al.ll201ll ). The survey samples a total volume of 
~1 Gpc^ over 816 deg^ of sky across seven fields, with 
an average target density of 350 galaxies per deg^. Pho- 
tometric data used in target selection for the WiggleZ sur- 
vey are f rom the the Galaxy Evolution Explorer (GALEX) 
satellite (jMartin et al.ll2005l ). the fourth data release of the 
Sloan Digital Sky Survey fSDSS: lAdelman-McCarthv et al] 
12006 ) in the northern galactic region, and the Canada- 
France-Haw aii Telescope S e cond Red-sequence Cluster Sur- 
vey (RCS2; iGilbank et al.l l|20lj )') in the southern galactic 
region. 

Whilst the WiggleZ galaxy redshift distribution peaks 
at 2: ~ 0.6 a long thin tail of high-redshift galaxies reaches to 
z — 1.5. These high-redshift galax;ies are identified by strong 
[Oil] emission lines, with the [Oil] doublet marginally re- 
solved, and have star formation rates of > 5OM0 yr~^. 

Galaxies selected for this study were drawn from the 
WiggleZ survey subject to a variety of criteria. Primary se- 
lection required (1) redshifts of 1.2 < z < 1.5, where the [O 
II] emission line falls between OH sky residuals at ~8000- 
10000 A, and (2) [Oil] emission exceeding 50 x 10" ^'^ ergs 
s~^ cm~^. We found ~300 targets in the WiggleZ survey 
that satisfied the aforementioned criteria. To reduce contam- 
ination from active galactic nuclei (AGN), starbursts were 
preferentially selected by (3) omitting galaxies with [Mgll] 
emission visible in the Wig gleZ spectrum which may be in- 
dicative of AGN activity (JGrandi fc Phillips! 1 1979? ) . With 
the additional criteria the number of targets was reduced 
to ~100 galaxies. Finally, for IFU observations with laser 
guide star adaptive optics (LGSAO), a suitable tip-tilt star 



is required to correct aberrations in the wavefront. We re- 
quire a star with apparent r-band magnitude rrii < 17 with 
an angular separation < 60 arcsec from the science target 
for optimal acquisition. This last criteria reduced our final 
target sample to ~50 targets. 

2.2 OSIRIS Observations 

Out of the possible targets, spanning across seven fields, 
we observed 17 star-forming galax;ies with OSIRIS (OH 



Suppressing InfraRed Imaging Spectrogr aph: [ Larkin et al 



2006i) with the Keck II LG SAO system l|Wizinowich et al 



l2006l : Ivan Dam et al.ll2006h . These 17 targets were prior- 
itized for having the brightest tip-tilt stars with shortest 
separation, and for which [OH] and Ha fall fortuitously in 
clean sky regions. OSIRIS is a lenslet array spectrograph 
with a 2048 x 2048 Hawaii-2 detector and spectral resolu- 
tion R ~ 3600. Observations were performed in September 
2008, May, September, and December 2009, and March 2010. 
The filters Hnl, Hn2, and Hn3 were used with central wave- 
lengths of 1.500, 1.569, 1.635 /^m respectively in the 0.05" 
pixel scale. The FWHM of the tip-tilt stars ranged from one 
and a half to two pixels, 0.062-0.1", with an average Strehl 
of 30% estimated from the tip-tilt stars. Seeing ranged from 
0.4" to 1.4". The galax;ies observed on these nights are listed 
in Table [T] 

The standard observing procedure was as follows. We 
first acquired the tip-tilt star and applied the optimal po- 
sition angle of OSIRIS to position the star within the un- 
vignetted field-of-view of the LGSAO system. Short, 60s, 
integrations were taken on the star for calculations of the 
PSF and to centre the star in the field. A blind offset was 
then applied to the telescope to acquire the target. The Wig- 
gleZ survey selected targets from S DSS and RCS2 which fo l- 
low different astrometry systems (jDrinkwater et al.ll201Cl ). 
In calcul ating offsets from the reference star t o the target, 
2MA SS (|Skrutskie et al.ll2006l ) and USNO-B (JMonet et all 
120031 1 coordinates were used for SDSS and RCS2 galaxies 
respectively. The typical science observing sequence was to 
expose on target for 900s in each of two nod positions, A 
and B, separated by 1.6" on the detector. If Hq emission 
was detected in the difference of the first two exposures, we 
continued to expose on the galaxy for --^1 — 2 hrs as neces- 
sary to achieve a high signal-to-noise detection. Multiple A 
and B positions were dithered by 0.05" around the base po- 
sitions in each exposure to remove bad pixels and cosmic ray 
contamination. Of the 17 observed targets. Ha was detected 
in the spectra of 15 targets. The two non-detections were a 
result of the [OIIIJA5007 emission line being mis-identified 
as [OII]A3727 in the WiggleZ spectrum. Of the 15 detected 
targets, observations for two galaxies are unusable for sci- 
entific reduction due to poor AO correction as a result of 
stray light on the low bandwidth wavefront sensor. The Ha 
luminosities of the final 13 galaxies are shown in Figure [T] 
with [OH] luminosities of all possible targets indicated by 
the small gray filled circles. 

Data reduction was completed with the OSIRIS data 
reduction pipeline versions 2.2 and 2.3 and custom IDL rou- 
tines developed for faint emission-line spectra. The pipeline 
removes crosstalk, detector glitches, and cosmic rays be- 
fore it mosaics individual exposures and assembles a re- 
duced data cube with two spatial and one wavelength di- 
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Figure 2. OSIRIS integrated spectra of each galaxy over all pixels satisfying the criteria outlined in Section 2. The vertical dotted lines 
denote the position of the redshifted locations of nebular emission lines: Ha, [NIIJAA6548, 6584, and [SIIJAA6716, 6731 at the systematic 
redshift. The solid red line is the best-fitting single Gaussian fits to each detected emission line. 
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Figure 3. OSIRIS Ha kinematics. Left: Ha flux map, middle: velocity map from the systemic velocity, right: velocity dispersion map. 
The black crosses mark the peak(s) in Ha flux. Velocity shifts from the systemic redshift and velocity dispersion are measured in km 
s~^ . Compass at the bottom right of the left panels show the true orientation. Letters at the bottom left of the left panels represent 
morphology classifications; M=Multiple emission, E=Extended emission, S=Single emission. Objects are presented in the order based 
on the morphological classification described in Section 3.1. 
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Figure 3. cont. 
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Table 2. Emission Line Properties 



-fh 

■'[Nil] 



WiggleZ ID 



ID 



/hq 



(10**^ ergs s-i) 



12 + log(0/H) 



R03 J032450240- 13550943 
S15J145355248-00320351 
R01J005822757-03034040 
R03J032206214-15443471 
S15J144102444-I-05480354 
R00J232217805-05473356 
S09J091517481-I-00033557 
S00J233338383-01040629 
S15J142538641-I-00483135 
S11J101757445-00244002 
S09J090933680-I-01074587 
S11J110405504+01185565 
S09J090312056-00273273 



WK0912_13R 


1.2873 


48 ± 16.1 


<8.4 


46.8 ± 13.2 


<8.47 


WK1002.61S 


1.3039 


81.5 ± 5.3 


<9.4 


81.9 ± 13.7 


<8.37 


WK0909.02R 


1.3616 


41.5 ± 7.9 


9.3 ± 4.4 


46.4 ± 15.2 


8.53 ± 0.29 


WK0909.06R 


1.4602 


94.9 ± 6.7 


22.8 ± 3.7 


126.2 ± 18.1 


8.55 ± 0.1 


WK0905.22S 


1.2810 


105.4 ± 4.3 


10.8 ± 2.3 


101.4 ± 13.1 


8.34 ± 0.12 


WK0912.01R 


1.2960 


57.3 ± 7.5 


<6.4 


56.7 ± 13.5 


<8.36 


WK1002.41S 


1.3330 


64.9 ± 4.3 


13.7 ± 2.4 


68.9 ± 14.4 


8.52 ± 0.1 


WK0809.02S 


1.4521 


74.7 ± 13.6 


6.3 ± 7.3 


98.0 ± 17.8 


8.29 ± 0.68 


WK1002.18S 


1.3044 


58.0 ± 4.8 


<7.6 


58.4 ± 13.7 


<8.4 


WK1002.14S 


1.3048 


92.8 ± 4.3 


20.6 ± 2.4 


93.5 ± 13.7 


8.53 ± 0.07 


WK1003.46S 


1.3074 


50.0 ± 4.5 


<5.6 


50.6 ± 13.8 


<8.36 


WK1003.52S 


1.3191 


86.1 ± 6.0 


20.4 ± 3.4 


89.0 ± 14.1 


8.54 ± 0.1 


WK0912_16S 


1.3260 


38.6 ± 4.4 


7.5 ± 2.4 


40.5 ± 14.3 


8.49 ± 0.19 



° Systematic redshift measured from Gaussian fit to Ha emission in integrated spectrum. 
'Emission line flux in units of 10"^'^ ergs s~^ cm~^. 



mension. First-order sky subtraction was achieved by the 
spatial ABBA nodding on the sky. Further sky subtraction 
was apphed using custo m IDL routines which employ the sky 
subtraction methods of IPaviea (|2007|). Galaxies were flux 
calibrated using telluric standard stars observed throughout 
the night. On average three telluric standards were observed 
each night. 



2.3 Kinematic Mapping 

In order to create 2-D spatial and kinematic maps of the Hq 
emission line, the sky-subtracted data cube was convolved 
with a Gaussian kernel with FWHM — 150 mas (200 mas 
in cases of low signal) to achieve a higher signal-to-noise ra- 
tio in each spaxel. A Gaussian profile was then fitted to the 
Ha emission in each spaxel of the cube with the wavelength, 
intensity and width as free parameters. Errors were calcu- 
lated by adding in quadrature an estimation of the sky noise 
from the variance of the spectrum offset from the emission 
lines with the Poisson error of the photon counting statis- 
tics in the emission lines. Spaxels were masked from the Ha 



Table 3. IRIS2 Observations of WiggleZ Galaxies 



ID 


Obs. 


J Exp. 


K Exp. 


J"^ 


K" 




Date 


timc(s) 


time(s) 






WK0912.13R 


2010/01 


1350 


4050 


20.71 


20.13 


WK0909.06R 


2010/01 


2700 


3360 


20.77 


19.98 


WK0905.22S 


2009/05 


1350 


1350 


21.36 


20.52 


WK1002.41S 


2010/01 


1350 


3600 


21.30 


20.57 


WK0809.02S 


2010/07 


1950 




20.66 




WK1002_18S 


2009/05 


1350 


2700 


20.84 


21.02 


WK1002.14S 


2009/05 


1350 


1350 


20.69 


20.20 


WK1002.52S 


2010/01 


1350 


3300 


19.94 


18.92 


WK0912.16S 


2010/01 


1260 


4050 


21.79 


21.93 



"Apparent AB magnitudes. The average error on J and K 
magnitudes are 0.11 



emission map and kinematic maps if the peak of the fitted 
Gaussian was less than one standard deviation above the 
median of the spectrum (corresponding to minimum signal- 
to-noise of ~2), or the width of the Gaussian was less than 
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the resolving power of the instrument. For spaxels surviving 
these cuts a flux was calculated from the fitted Gaussian. 
The Ha flux images represent this array of fluxes. The inte- 
grated spectra, shown in Figure [2] were created by summing 
all unmasked spaxels in the spatial directions. A Gaussian 
was fitted to the observed-frame integrated spectrum to ob- 
tain the total Ha flux and the systematic redshift, Zays given 
in Table (2] A velocity map was constructed for each galaxy 
using the velocity shifts of the emission line in each spaxel 
relative to Zsys . The velocity dispersion was measured in each 
spaxel from the r.m.s. of the fitted Gaussian line profile cor- 
rected for instrumental broadening. All kinematic maps are 
shown in Figure [S] 



2.4 NIRSPEC Observations & Data Reduction 

When the Keck LGSAO system was unavailable due to bad 
weather, technical faults or space events, observations of 
available targets were tak en with the near-inf rared echelle 
spectrograph (NIRSPEC: iMcLean et al.lll99^ to measure 
the Hj3 and [OHI] emission line ratio. The ratios of emission 
lines from these observations, obtained for VK-ft'0905_22S, 
WKim2_19>S, WKim2J,2S, WKim23lS, are used to 
constrain their AGN contamination as discussed in Section 
3.4. 

The NIRSPEC detector is a 1024x1024 ALADDIN-3 
InSb detector with spectral resolution R ~ 2000. All obser- 
vations were carried out in February 2010 with the 4-pixel 
slit (42x0.76") in the NIRSPEC-2 filter (f.089-1.293 /im). 
Seeing in K was 0.8" on average. 

The typical observing sequence was to acquire a guide 
star, typically the star already selected for use with OSIRIS 
LGSAO system. Short (~I0s) images were taken with the 
NIRSPEC guider to center the star at the central position 
on the slit. Once centered a blind offset was applied to the 
star position. Longer images were taken (2 x 20s) at the tar- 
get position. If the target was observed in these images it 
was centered on the slit by nudging the slit left and right 
in single pixel intervals. Once centered spectroscopic object 
exposures were taken with typical sequence of 2 x 600s in 
each of two nod positions, A and B in the sequence ABBA. 
Occasional shifts of a few pixels were needed to keep the 
object on the slit throughout the object exposures. No con- 
tinuum was detected and the ratio of [OIII]A5007/H/?A486f 
was calculated from the raw data. 



2.5 IRIS2 Observations 

Near-infrared imaging was obtained for 9 out of 13 galaxies 
in the s ample with the Inf rared Imager and Spectrograph 
(IRIS2, iTinnev et al] 120041 ') on the Anglo- Australian Tele- 
scope (AAT). J and Ks bands, with central wavelengths 
Ac = 1.245 /im and Ac = 2.1444 /im respectively, were chosen 
in order to image the older stellar population of each galaxy 
to obtain robust estimates of stellar masses. The IRIS2 de- 
tector is a 1024 x 1024 Rockwell HAWAILl HgCdTe infrared 
detector, with a scale of 0.4487 arcsec pixeP^ and field of 
view of 7.7 arcmin^. Observations were carried out during 
May 2009, January 2010, and July 2010. Seeing in Ks ranged 
between 0.9"-2.5", 0.9"-1.6", and 1.3"-2.0", respectively. 
The typical J exposure sequence was a series of 15 sec- 



ond exposures totaling 1350 seconds on the target. The typ- 
ical Ks exposure sequence was a series of 10 seconds expo- 
sures totaling 1350 seconds on target. Most galaxies required 
2 or 3 sequences in Ks in order to obtain a detection with a 
signal-to-noise ratio greater than 5. Due to the time-sensitive 
variation of the infrared night sky, we used a dithering strat- 
egy of successive observations with 9 dithers per sequence 
with a maximum dither length of 20 arcsec. 

Data reduction and mosai cking was carried out b y data 
reduction pipeline ORAC-DR (|Cavanagh et al.l [20031 '). Each 
image was flattened with the other images closest in time 
and position in the dither pattern. The dithering script av- 
eraged out bad pixels on the detector and applied a nearest 
neighbor smoothing. 

Magnitudes were calculated in the AB system using 
aperture photometry with a circular aperture of 4" with 
a concentric sky annulus. A zero point was derived for each 
reduced mosaicked image by comparing the measured aper- 
ture counts with known 2MASS magnitudes for several stars 
in the image. The zero point was then used to photometri- 
cally calibrate the image. On average, four uncrowded cen- 
tral stars with ii'-band magnitudes of ii' ~ 15 were selected 
to derive the zero point. If the standard deviation of the zero 
points from these stars was greater than 0.3, the stars pro- 
ducing the most deviant zero points were replaced with other 
stars in the fleld. Errors were calculated from the variance in 
the sky aperture. The sky background was calculated using a 
clipped mean. The IRIS2 system most closely resembles the 
Mauna Kea Observatory (MKO) system. Equations [T] and 
[5] were used to transform between the MKO and 2MASS 
systems. 



Jmko ~ 


■- -0.03 - 


-0.03 X (J- 


^Ks) + J 


(1) 


Kmko ~ 


: -0.01 - 


-0.01 X (J- 


-Ks) + K 


(2) 



where J and Ks are given in the 2MASS system and 
Jmko and Kmko are given in the Mauna Kea Observatory 
system used by IRIS4J. Derived values of J and ii'-band 
magnitudes are given in Table [S] 



2.6 Stellar Mass Fitting 

Stellar masses and star formation rates were calcu- 
lated using an spectral energy dis tribution (SED) fit- 
ting routine (jGlazebrook et al.l 120041 ') which employs the 
Baldry fc Glazebrook initi al mass function (IMF) (BG03; 
iBaldrv fc Glazebrookl [200^ ) . The BG03 IMF fits local cos- 
mic luminosity densities well, with a high-mass slope similar 
to the clas sical Salpeter IMF but with a more realistic break 
at O.5M0. iHopkins fc BeacomI (|2006r ) have shown that the 
BG03 IMF gives a better agreement between cosmic star- 
formation histories and mass assembly. 

Infrared JK data from IRIS2 was combined with multi- 
band photometry available in WiggleZ flelds, including ugriz 
from SDSS, grz from the RCS2 Survey, and NUV from 
GALEX. Tw o-component modeling was adopte d (using PE- 
GASE.2 from lFioc fc Rocca-Volmerang3 (|l997l ') to calculate 
spectra) to assess the possible biases from starbursts on the 
calculated masses and to allow a range of dust extinction (0 



http://www.aao.gov.au/iris2/ 
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Table 4. Ha Kinematic & Morphological Properties 



ID 


shear 


"mean 


<ct 


^shcar/ '^mcan 


'^1/2 


Class<= 




(km s-i) 


(km s-i) 


(kms-i) 




(kpc) 




WK0912_13R 


130.5 ± 11.0 


81.6 ± 27.7 


104.3 ± 27.1 


1.6 ± 0.6 


3.8 


M 


WK1002_61S 


71.5 ± 10.7 


85.1 ± 20.3 


88.4 ± 4.4 


0.9 ± 0.2 


3.8 


M 


WK0909.02R 


136.4 ± 6.1 


87.9 ± 21.9 


103.7 ± 15.5 


1 ± 0.3 


2.6 


M 


WK0909.06R 


160.8 ± 2.4 


92.6 ± 27.8 


153.3 ± 8.4 


1.6 ± 0.5 


4.2 


M 


WK0905_22S 


55.5 ± 6.5 


98.2 ± 10.9 


107.4 ± 3.3 


0.6 ± 0.1 


2.3 


E 


WK0912_01R 


99.4 ± 7.1 


98.6 ± 19.6 


112.1 ± 11.0 


1.3 ± 0.3 


2.8 


E 


WK1002.41S 


80.2 ± 5.4 


91.4 ± 13.0 


116.3 ± 5.9 


1.7 ± 0.3 


3.1 


E 


WK0809.02S 


48.5 ± 9.4 


92.7 ± 16.6 


97.6 ± 13.8 


0.4 ± 0.1 


2.1 


E 


WK1002_18S 


81.1 ± 17.5 


102.1 ± 20.0 


106.6 ± 6.8 


0.7 ± 0.2 


2.1 


S 


WK1002_14S 


48.4 ± 9.5 


113.7 ± 8.8 


133.3 ± 4.8 


0.5 ± 0.09 


2.1 


S 


WK1003.46S 


65.6 ± 16.7 


124.8 ± 24.9 


138.6 ± 9.6 


0.5 ± 0.2 


1.8 


S 


WK1003.52S 


73.7 ± 18.2 


125.4 ± 20.9 


145.4 ± 8.1 


0.7 ± 0.2 


1.8 


S 


WK0912.16S 


54.9 ± 12.9 


85.1 ± 18.6 


104.4 ± 9.1 


0.7 ± 0.2 


2.5 


S 



"Velocity shear measured across the Ha velocity map as described in Section 3.2. 

'Flux-weighted mean velocity dispersion calculated from a Gaussian fit to the spectra of individual spaxels of the OSIRIS Ha detection. 

'^Velocity dispersion calculated from a Gaussian fit to the integrated spectrum of the OSIRIS Ha detection. 

''Half-light radius of OSIRIS Ha emission in circular aperture. 

'^Morphological classification; EiExtended emission, MiMultiple emission, S:Single emission. 



^ Av sC 2 mag) and metallicity (0.0004 ^ Z ^ 0.02). The 
primary component was modeled using SFR oc exp(-t/r) 
with r =0.1, 0.2, 0.5, 1, 2, 4, 8, & 500 Gyr (the first approx- 
imates an instantaneous starburst and the last a constant 
SFR). The full distribution function of allowed masses was 
calculated by Monte-Carlo re-sampling of the photometric 
errors. The final masses and error bars, given in Table O 
represent the mean and standard deviation of this full dis- 
tribution function. In the cases where observed UV, optical 
and NIR magnitudes were available the errors in the de- 
rived stellar masses are on average 5[log(M, [M©])] — 0.15. 
When NIR magnitudes were not available, the errors average 
<5[log(M4M0])]= 0.62. 

We do not us e the models of lMaraston et al.l l|2006l ) as 
iKriek et all (|2010t ) show a low contribution from TP-AGB 
stars for post-starburst galaxies a t high redshift. For mor e 
information on the SED fitting see lClazebrook et al.l (|200J). 



3 RESULTS 

In Figure [3] we present Ha spatial (left panels) and kine- 
matic maps of our 13 galaxies observed with OSIRIS with 
spatial extent ranging from 2 — 12 kpc. A broad range of 
morphologies and kinematics are observed with 8/13 galax- 
ies showing some velocity shear and the remaining showing 
disturbed kinematic morphologies (middle panels). The dis- 
persion maps for the entire sample are chaotic with peaks 
in velocity dispersion exceeding 100 km s~^ (right panels). 
In the next subsections we will discuss the Ha morphology 
and kinematics, emission-line properties, and masses of the 
sample. For convenience IDs are shortened and referred to 
by their last three unique digits. 

3.1 Classification 

We initially classified the galaxies into three groups based 
solely on Ha morphology: single emission, extended emis- 
sion, and multiple emission. Galaxies classified as single- 
emission have a single resolved source of Ha emission: 16S, 



14s, 18S, 46S, 52S. Systems classified as extended emis- 
sion show multiple unresolved or extended components: 02S, 
22s, OIR, 41s. Systems 02R, 06R, 13R, 61S are classified 
as multiple-emission galaxies, which have multiple, resolved 
regions of Ha emission. The centres of the separate compo- 
nents are marked in Figure [3] These classifications are made 
for the purpose of discussing properties of galaxies with sim- 
ilar morphologies. The classifications are used to order Ta- 
bles and Figures. Kinematic classifications are discussed in 
Section 4. 



3.2 Kinematic Quantities 

A flux-weighted velocity dispersion (crmoan) was calculated 
from the dispersion and flux across all spaxels. We define 
velocity dispersion by amean, which is representativ e of the 
inter nal motions of the gas (jLaw et al. 2007 ; Green et al.l 
l2010f ). For comparison, a line-of-sight velocity dispersion 
(o"nct) was calculated from the width of a Gaussian fit to 
the Ha emission of the integrated spectrum. The measured 
values of (TniGan are consistent with or slightly less than cTnct, 
as expected from velocity smearing. The values of amcan in 
this sample reach to ~ 130 km s~^ , and are ty pically higher 
than in the samples of iLaw et al.l (1200 !/) and Epinat et al.l 
(2009). We note that the 0.15" smoothing described in Sec- 
tion 2.3 produces a -1-10 km s~^ median boost compared to 
Cmoan Calculated from the non-smoothed data which is less 
than the statistical error. We give the (Jmoan values from the 
smoothed data in Table |4] for comparison to the literature. 
We do not make a correction based on beam smear- 
ing resultant from the PSF, which can be a problem for 
interpreting kinematic observations at high redshift (e.g. 
Law et al.ll2009l:IWright et alllJOOgl : iForster Schreiber etal] 



2OO9I : lEpinat et al.l |2010|) . Although all observations were 
taken in LGSAO mo de, with a PSF of FWHM~850 pc, 
iLaw et all (|2009l ) and IWright et al] (|2009t ) note that beam 
smearing can still have the effect of reducing the maximum 
rotation velocity and infiating the velocity dispersion for 
edge on disks in particular. However, as the inclination is 
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Table 5. Stellar Properties 



ID 


log[M.(M0)] 


Av 


SFRh^ 


SFRfoii] 


SFRsed 


SSFR 






(mag) 


(M© yr-i) 


(Mo yr-l) 


(M© yr-i) 


(10-9yr-i) 


WK0912_13R 


10.7 ± 0.2 


1.1 ± 0.3 


20.3 ± 6.8 


43.3 ± 0.9 


116 ± 40.4 


1.8 


WK1002_61S 


10.7 ± 0.7 


0.2 ± 0.2 


35.6 ± 2.3 


69.0 ± 2.8 


25.6 ± 13 


2.9 


WK0909.02R 


10.3 ± 0.6 


0.7 ±0.4 


20.2 ± 3.8 


33.1 ± 0.8 


57.4 ± 41.5 


4.6 


WK0909.06R 


11.0 ± 0.2 


0.8 ± 0.3 


54.8 ± 3.9 


120.3 ± 14.9 


91.8 ± 48.7 


2.09 


WK0905_22S 


10.3 ± 0.2 


0.3 ± 0.2 


44.0 ± 1.8 


202.4 ± 23.8 


44.6 ± 15.1 


10 


WK0912_01R 


10.0 ± 0.6 


0.5 ± 0.2 


24.6 ± 3.2 


59.0 ± 3.8 


40.3 ± 13 


11.2 


WK1002.41S 


10.2 ± 0.1 


0.4 ± 0.1 


29.9 ± 2.0 


112.8 ± 5.7 


60.3 ± 9.3 


8.1 


WK0809.02S 


11.1 ± 0.6 


0.0 ± 0.1 


42.5 ± 7.7 


29.7 ± 0.8 


19.5 ± 6.9 


1.5 


WK1002_18S 


10.2 ±0.1 


0.5 ± 0.1 


25.3 ± 2.1 


188.6 ± 14.2 


64 ± 17.4 


6.09 


WK1002.14S 


10.6 ± 0.2 


0.3 ±0.1 


40.6 ± 1.9 


43.9 ± 0.7 


55.5 ± 16 


4.5 


WK1003.46S 


9.8 ± 0.6 


0.2 ± 0.3 


22.0 ± 2.0 


41.3 ± 4.8 


22.4 ± 17.1 


12.8 


WK1003.52S 


11.7 ± 0.1 


0.0 ± 0.1 


38.6 ± 2.7 


103.7 ± 16.9 


15.9 ± 2.2 


0.3 


WK0912.16S 


9.8 ± 0.2 


0.1 ± 0.1 


17.6 ± 2.0 


64.1 ± 6.8 


12.6 ± 3.2 


12 



''Errors on the [Oil] SFRs are statistical errors only. 
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Figure 4. Flux weighted v elocity dispersion as a function of 
L[jc from lGreen et al.l 1 120101) with galaxies from this sample as 
red diamonds. The WiggleZ points extend this relationship to 
higher o-mcan and Ha luminosity. 



not well constrained in this study it is unclear when and 
what correction would be appropriate. The same authors 
note that beam smearing does not explain the high veloc- 
ity dispersions of objects which show no rotational velocity 
structure. 

The observed dispersion and luminosity of the WiggleZ 
sample conti nues the rela ti onshi p, spanning over < z < 3, 
observed by I Green et al.l ([20101) for star-forming galaxies, 
recreated in Figure |4l The galaxies in the WiggleZ sample 
follow this relationship tightly with the highest dispersions 
and H g luminositie s of th e IFS samples introduced in Fig- 
ure [l] iGreen et al.l (120101 ) argue that the correlation seen 
implies star formation rate is the driving factor of the high 
velocity dispersion, not mass or gas fractions. Consequently, 
arguments that high velocity dispersions at high redshift are 
evidence for cold-flow accretion are weakened as the cold 



Figure 5. Logarithmic velocity shear as a function of logarith- 
mic st e llar mass for our 13 galax i es and galaxies fromlLaw et al.l 
ll2009l):[Gen2el et a.1.1 l|2006l , l2008l) : lForster Schreiber et al.1 1|2006| ) 
and lEpinat et al.l (|2009|). Data are not corrected for inclination. 
All stellar masses are converted to the Salpeter IMF for ease of 
comparison with the literature. A median error for the WiggleZ 
points is plotted in red in the lower right corner. The red dotted 
line is the best-fitting linear model to the 13 galaxies in our sample 
[log(t;ghear /km s-i) = 0.03 X log(M*/MQ)+1.58]. The blue dot- 
dash line is the best-fitting linear model to the data from all the 
samples [log(Dshcar /km s~l) = 0.25 X log(Mt/M(7)) -0.84], the 
blue d otted line is the best-fitting linear model from iLaw et al.l 
&i) [log(«shear /km s"!) = 0.29 X log(M. /Mq )-1.25] and the 
black dashed line is at the value of (log{crnioan)). 



flows are no longer in place at the current epoch to produce 
the high dispersions observed in local galaxies. 

A velocity shear was mea sured acro s s eac h galaxy. 



^shca 



= i(«„ 



), following iLaw et al.l (|2009l ). In order 



to prevent biasing from outliers, iimax and Vmin are calculated 
from the mean of th e highest and lowest 5 % of the veloc- 
ity map respectively (|Goncalves et al.ll2010l ). The values of 
^shoar Underestimate the true rotational velocity of the gas 
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in the cases that the galaxy extends to greater radii than 
detected by OSIRIS. The values of amcan, crnct, Wshcar, and 
f shear /cTmcan are given in Table [H We do not find a corre- 
lation between velocity shear and stellar mass in this sam- 
pl e of 13 galaxies, in disagreement with the z ~ 3 sample 
of iLaw at al.l (120091) an d Lyman break analogs at z ~ 0.1 
JGoncalves et al.l 1201(3 ). However, in the context of other 
star-forming galaxies at 1 < 2 < 3 a positive correlation is 
seen, as shown in Figure \S\ The interpretation of the rela- 
tionship is that stable rotation is more prevalent in galaxies 
that have assembled a relatively large stellar population. 
However, the addition of this sam ple does n o t affe ct the 
correlation originally published by iLaw et al.l (|2009|) . It is 
unclear if the velocity shear from all samples are compara- 
ble given the different spatial scales probed by OSIRIS and 
SINFONI and the different signal-to-noise reached for each 
set of observations. 



3.3 Star Formation Rate Indicators 

Star formation rates were measured from WiggleZ [Oil] and 
OSIRIS Hq luminosities. Due to uncertainties arising in the 
flux calibration of optical fiber spectroscopy, the WiggleZ 
[OH] fluxes are only accurate to within 40% due to system- 
atic errors. OSIRIS Ha flux measurements are comparably 
more accurate and are thus used as the primary measure of 
the nebular star formation. Star f ormation rates (SFR) were 
calculated from iKennicutj (|l998l ). 



SFRh<.[Mo yr- 



SFRoii[M0 yr" 



1 Ljia [ergs s ^] 
1.82 1.26 X 10*1 
1 Loii [ergs s~^] 



(3) 
(4) 



1.82 7.1 X 10*0 

where the 1/1.82 converts to the BG03 IMF. 

Further measurement of star formation is derived from 
the SEDs generated from the stellar mass fitting described 
in Section 2.5 (SFRsed). A dust extinction, Avseo, is es- 
timated from the SED fitting described in Section 2.6 and 
given in Table [S] The mean error on Avsed is 0.2. 

3.4 Metallicity 

[Nil] was detected in 8 of the 13 galaxies in our sample. 
The detections and upper limits are used to investigate the 
metallicity and level of contamination from hosting an AGN 
by using the [Nil] A6584/Ha rati o. Metallicity was calculated 
from this ratio, as calibrated bv lPettini fc Page! (|2004l ). 



12 + log(g) 



8.90 + 0.57 xlog(l^) 



(5) 



^H' °' Ha 

and given in Table O The metallicities of the sample are 
all sub-solar and are co nsistent wit h the results of IFS 
sampl es at z ~ 1 — 3 of iGenzel e t_ al. (2008); Stark et al.l 
ll2008l): Haw et al.1 (|2009l ): iForster Schreiber et all (|2009l ): 
IWright et al.1 (|2009|) and DE EP2 star-forming galaxies at 
z ~ 1.3 (|Shaplevet ai]|2005l ). The typical errors in 12 + 
log(§) are 0.2. 

Metallicity measurements calculated from Equation [5] 
alone can be biased due to elevated levels of [Nil] caused 
by AGN contamination. We look for this bias by con- 
straining the position of the galaxies in this sam- 
ple on the emission line diagnostic, or BPT, diagram 
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Figure 6. The emission lin e diagnostic diagram of 

iBaldwin. Phillips fc TerlevichI lll98ll ). The black points are Wig- 
gleZ galaxies and the red diamonds are 52S, 22S, 61S and 18S 
from this sample. The shaded gray area is the range contain- 
ing the [Nil] /Ha ratio for the remainder of the systems in 
this sampl e. The dashe d curved line is the empirical line from 
JKauffman n et al.l 1120031). the so lid curved line is the theoretical 
curve from iKewlev et al.l | |200J) and the vertical and horizontal 
lines are used to identify Seyferts and LINERS as labelled. The 
WiggleZ sub-sample presented here is primarily in the region of 
the diagram were star formation is the dominant mechanism driv- 
ing line emission. 



(JBaldwin. Phillips fc TerlevichI [l98ll ), shown in Figure H 
NIRSPEC observations of [OIII]A5007, 4959 and H/3 of four 
galaxies in the sample further quantify AGN contamination 
using the [OIII]A5007/H^ ratio. Two of the [Nil] /Ha mea- 
surements are upper limits (see column 5 of Table 2), and 
in all cases the [Nil] /Ha ratios are significantly lower than 
the level expected for the given values of [OIII]A5007 / H/3 
if there was a significant component of AGN emission. The 
remainder of galaxies all lie to the left of log([NII]/Ha) = 
—0.6, in the shaded gray region of Figure (6] Although AGN 
contamination can not be ruled out completely for these 
sources, their integrated nebular emission is clearly domi- 
nated by star formation. 



3.5 Mass Estimates 

The sample stellar masses range from 9.8 < log(M*[M0]) < 
11.7, comparable to masses from other IFS s amples cover- 



ing the same redshift range (1.2 < z < 1.7; Epinat et al.l 



I2OO9I : I Wright et aDl2009l : iForster Schreiber et al.ll2009t ) but 



with a higher (M«) than high-reds hift IFS samples [z > 2; 
iLaw et al.ll2009l : lGenzel et al.ll2008l ). No correlation is found 
between Ha morphology and mass; the three most massive 
systems cover all three morphological classes, and the three 
lowest-mass galaxies cover two morphological classes. 

With stellar mass estimates we are able to identify 
where these systems lie relative to the main sequence of 
star formation (SFR vs. M*) as a function of redshift 
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Figure 7. Ha star formation rates of IFU samples at 2; > 1 as a function of stellar mass. Stellar masses and SFRs for all samples are 
converted to the Salpeter IMF for ease of comparison to the literatur e. The symbols ar e as described in Figure [T] The solid lines are 
from the Great Observatories O rigins Deep Survey f GOODS) given in J Elbaz et al.l 1 1200711 and lDaddi et al.l 1(2003) • The dotted lines are 
from the Millennium simulation dSpringel et al.ll2005l : ICroton et al.ll2006l ) at redshifts 2 = (light blue) and 2 = 1 (blue). An extinction 
correction is made from E(B— V)sbd- 



informing on the i ntensity of their star formation (e.g. 
iNoeske et al] l2007l : iRodighiero et all I2OI0I '). If positioned 
above the main sequence, outside the hmited range of SFRs 
at a given M,, the star formation of the sample may be 
indicative of episodic variations rather than being associ- 
ated with the larger population of continuous star-forming 
gala:xies. Furthermore, the main sequence evolves with red- 
shift in which hig h er SF Rs were more typical in past epochs. 
IRodighiero et al.l l|20ld ) show that the mean specific SFR 
(SSFR=SFR/M,) of gala^cies decreases by a factor of ~ 15 
from z = 2 to 2 = for the most massive galaxies. The main 
sequence for IFS samples is shown, as SFR vs M,, in Fig- 
ure [7] in two redshift bins to account for this evolution. The 
SSFR is calculated, using SFRhq, and presented in Table (5] 

In the lower-redshift bin of Figure[71 1.2 < z < 1.7, the 
star- fo rming gala:xies of the SINS survey and IWright et al.l 
([2003) form a sequence of star formation which scatter be- 
tween the GOODS z = \ and z = 2 empirical relations. The 
WiggleZ points in contrast do not show an increase in SFR 
with mass and therefore the galaxies in this sample below 
M* = 5 X 10^° M0 scatter above the main sequence and the 
galaxies with masses > 5 x 10^" M© scatter below the main 
sequence. We note that the extinction correction applied 
here is derived from the SED fitting and future measure- 
ments of the dust content from far-infrared indicators will 
provide a more accurate correction affecting the location of 
this sample along the main sequence. 

The high-redshift bin, 2.0 < z < 2.7, shows a sim- 
ilar main sequen ce from the star- forming galaxies of the 
SINS Survey and iLaw et al.l ((20091) which scatters around 



the GOODS z = 2 empirical relation. The mean SFR in the 
low-redshift bin ((z) ~ 1.52) is a factor of 0.7 less than 
the the mean SFR in the high-redshift bin ((2) ~ 2.3). 
Sub-millimeter galaxies (SMGs) are included in Figure [7] 
for comparison to the IFS samples. SMGs arc c onsidered 
as the 'maximum starbursts' (|Tacconi et al.ll2008l ) with in- 
frared determined SFRs > 1000 M© yr~^. For consistency, 
only SMGs with published nebular emission are included 
in the co mparison to IFS results (jSwinbank et al.l 12004 : 
Takata et al. 2006 ). SMG stellar masses are taken from 



Michalowski et al.l (|20ld ). SFRs for all samples are calcu- 
lated from measured nebular flux values with an extinction 
correction from Av,sed and with a Salpeter IMF for consis- 
tency. 

We find that the SMGs in the low-redshift bin are well 
matched statistically to the lower mass data presented here 
and higher on average than the other IFS samples. In the 
high-redshift bin the SMGs are offset in mass from the IFS 
samples. Unfortunately, the low number of SMGs with mea- 
sured line emission limits an in-depth comparison to this 
interesting galaxy population. However, initial comparisons 
to the available data indicate the WiggleZ sample presented 
here may be related to SMGs with more episodic star for- 
mation than LBGs at a given redshift and stellar mass. 

3.6 Disk Fitting 

The velocity maps of the galaxies were modeled to determine 
whether the data are consistent with disks. Inclined disk 
models are generated varying the kinematic centre (a;o,j/o), 
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Figure 8. Disk fitting results for tlie 7/13 galaxies well-fit by disks in the sample. The first panel shows the smoothed velocity map, 
the second panel shows the smoothed modeled velocity map fit to the unsmoothed data, the third panel shows the residual from the 
difference of panel 1 and panel 2 and the final panel shows the binned rotation curve across the major axis with the model overlaid in 
red. 
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Table 6. Disk Model Parameter Results 



"xT 



ID 



^0_14 

53±2 
287ti° 
226121 

219±14 
320±36 



Vp sin(i) 
(km s-i) 



^^ 



(kpc) 



l±l 
18«3 

1+6 



(km s^i) 



Residual 

(km s^i) 



WK0912. 
WK0909. 
WK0905. 
WK0912. 
WK1002. 
WK1002. 
WK1002. 



.13R 
.06R 

.22S 

.OIR 

.41S 

.14S" 

.18S" 



73 
187±7 
102 
83 
141 
42 
24 



+80 

26 

+49 

43 

126 

20 
+49 

30 
+32 

21 



10^ 



17 
16 



+3 
-17 
+8 
3 



8I9 


3.2 


1.97 


43l^ 


4.3 


5.12 


7+32 
-'-22 


3.5 


2.27 


14+7 


2.6 


3.44 


7+28 
"'-26 


1.5 


0.99 


I5I5 


1.1 


1.71 


28±7 


2.1 


2.58 



"Kinematic center fixed to peak of Ha flux distribution. 



the inclination (z, where i=0 is a face-on disk), the position 
angle (0), the systemic velocity (Vb), the plateau velocity 
(Vp), and the p lateau radius (rp) using the prescription of 
iBegemanI l| 19891) . 

V{x, y)^Vo + V,{r) sin(i) cos(6>), (6) 

where 14 is the circular velocity at a given radius from 
(a;o,J/o) and 9 is the azimuthal angle in the plane of the 
galaxy, defined as: 

-(a:: - a;o) sin(0) + {y - yo) cos{(l)) 



cos(6l) = 
sin(e) = 



-(x - 3:0) cos(<?!>) + {y - yo) sm((j}) 



r cos(i) 
The circular velocity is defined as: 



K(r) 



Vp 



r <r, 



Vp, r >ri 



(7) 



(8) 



for a flat rotation curve (|Wright et al.l 120071 '). The 'flat' 
disk model has been shown statistically to reproduce fit- 



ting parameters more successfully than the exponential disc, 
isothermal sphere , and arctangent m odels in artificially red- 
shifted IFS data |Epinat et al.ll201Cll l. The models are con- 
volved with the observed PSF, modeled from a tip-tilt star 
observation, and are fit to the data prior to the 0.15" 
smoothing to reduce the effects of beam smearing. The best- 
fitting models were determined by minimizing the reduced 
chi-squared value (xr) using Monte Carlo Markov Chain 
methods. The parameter results of the best-fitting models, 
Xr values, and the average residual velocity are given in 
Table |6l Example models are shown in Figure [S] The resid- 
ual velocity, Vias , is calculated from the pixel average of the 
absolute value of the difference of the data and model. 

The kinematic center is not well constrained in single 
emission galaxies, particularly in the single emission galax- 
ies that show no shear or have fewer spatial pixels. In these 
cases, specified in Table [l] the kinematic center is fixed to 
the peak of Ha emission. Although the peak of Ha emis- 
sion is not n e cessar ily the dynamical center of the galaxy, 
IWright et al.l (|2009l ) find that it is a reasonable approxima- 
tion. 
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We find that 7/13 of our galaxies {06R, OIR, 13R, 4IS, 
22s, 14s, 18S) are well fit by a disk model. To determine 
which objects are well fit by this simple disk model we con- 
sider the velocity residual, reduced chi-square value and a 
qualitative measure of coherent velocity shear. We note that 
Xr is biased as the spaxels are correlated due to the intrin- 
sic PSF and is not alone a reliable indicator. Furthermore, 
small fluctuations in the velocity map resulting from galactic 
winds and/or interactions between individual star-forming 
regions will inflate the residual and reduced chi-squared val- 
ues as they are not accounted for in the model. 

A known approximate degeneracy exists b etween V^ 
and i in the form Vp sin(t) in Equation |6] (e.g. iBegemanI 
1 19891 : iLaw et al.ll2009l : lEpinat et al.ll2009l ). In previous low 
signal-to-n oise IFS work the in clination i s set to a constant 
of i = 45 l|Wright et al.1 l2009t ) however, iLaw et al.1 (|2009l ') 
derive i — 57.3 as the correct average inclination between 
our line of sight to a given disk and the vector normal to 
the disk. The most robust method may be to fix the in- 
clination to values derived from deep broadband imaging 
IjEpinat et al.ll2009l ). however we consider VpSin(i) as one 
parameter in this study as we lack sufficiently deep imaging 
to determine it from photometry. 

We find another degeneracy present in 4^3, I4S and 
13R between the plateau radius and plateau velocity. This 
degeneracy arises from the linear nature of the velocity gra- 
dients observed and is likely a consequence of not probing to 
faint enough surface-brightness emission to detect the kine- 
matic turnover. A characteristic peak of velocity dispersion, 
a result of rapid unresolved change in velocity, would not be 
expected in the dispersion maps. Similarly, face-on and low 
mass galaxies which have a low velocity gradie nt show a very 
faint or no central peak of velocity dispersion (jEpinat et al.l 
[2OI0I) . As a result we do not use the central peak of velocity 
dispersion as an indicator for disk classification. 

As found in previous IFS studies of high-redshift galax- 
ies the Vshear/o"mGan (Table |4]) valucs of this sample are 
approximately unity. The v/a measurement is often used 
to make a split between rotationally supported systems 
{v/(T > 1) and systems dominated by dispersion (v/g < 1). 
In this sample 02R, 06R, OIR, 13R have v/g > 1 while 02S, 
22s, 16S, 46S, 52S have v/g < 1 and 4IS, 18S, 61S have 
v/g ~ 1 within the errors. This correlates with the results 
of disk modeling in which 06R, OIR, 13R and 4^3 are well 
fit by disks and 02S, 16S, 46S, 52S are not well fit by disks. 
However, this picture fails for 22S, I4S which are well fit 
by disks but have v/g < 1 and 02R which likely has an el- 
evated v/g as a result of an outflow rather than from the 
signature of rotational support (See Appendix). It is likely 
that the discrepancy between these two methods is a result 
of the detection limit of the observations. The Yghear mea- 
surement probes only the inner star-forming regions of the 
galaxies and is not representative of the overall rotational 
velocity of the galaxy, thereby reducing the value of v/g. 
Other factors are known to reduce the observed values of 
v/g such as t he foresh ortening due to the inclination by a 
factor of sini (jLaw et al . 2009). Because the inclination val- 
ues in this study are poorly constrained this correction is 
not applied. As a result, more weight is given to the results 
of disk modeling rather than v/g in classifying galaxies as 
being rotationally supported. 



3.7 Surface Brightness Profiles 

The effective radius was calculated by fitting convolved 1-D 
Sersic models to the Ha surface brightness profiles of single- 
emission galaxies. In order to fit model Sersic profiles, the 
central spaxel of the galaxy or star-forming complex was 
first identified from the peak of the fiux maps. The distance 
to the remaining spaxels was calculated from the central 
point assuming a circular profile, reducing the number of 
free parameters by setting the aspect ratio to unity and po- 
sition angle to zero. Individual measurements are averaged 
in radial bins. The models were convolved with the PSF. The 
best-fitting Sersic profile model was determined by minimiz- 
ing the chi-squared value (x^) with three free parameters, n, 
ro, and Eo, where n is the Sersic index, ro is the effective ra- 
dius and Eo is the surface brightness at the effective radius. 
We assume this simple model as the inclination and position 
angle determined in Section 3.6 are not well constrained. 

The fitted results are shown in Table [7] with retr rang- 
ing from 0.4 kpc < r^a < 1.9 kpc and n ranging from 
0.5 < n < 2.0. An example model fit to the data is shown 
in Figure [9] It is not completely understood what the fiux 
profile of Hq emission reveals about the parent galaxy as 
Sersic profiles are more commonly fit to broadband emis- 
sion representative of continuu m emission from continu- 
ous star formation (|Sersid 119681 ). However, the SINS sur- 
vey obtained deep rest-frame optical images of six irregular 
galaxies in their sample and found surprisingly good mor- 
phological agreement between Ha and broadband images 
(jForster Schreiber et al.ll201(J ). They find the Sersic param- 
eters of Ha emission and _H"-band emission (rest-frame op- 
tical) are consistent for the six galaxies with a difference on 
average of 0.08 for the Sersic index and agreement within 3% 
for reff. The Sersic indices are comparable to our Ha results 
with n < 1 in all cases however they find greater effective 
radii with reff > 2 kpc in all cases, most likely due to the 
coarser resolution from their non-AO corrected observations. 

The fitting method described above is applied to avail- 
able r-band images from RCS2 only, as galaxies are unre- 
solved in SDSS images. The RCS2 images, taken in aver- 
age seeing of 0.7", are sufficient to fit Sersic profiles to 3 of 
13 systems of the sample. We obtain results in the range 
0.2 < n < 1.0. A direct comparison to the Ha Sersic param- 
eters however is not possible as the only RCS2 data available 
is of extended or multiple emission galaxies. Deeper imaging 
of the whole sample is needed to reveal whether the under- 
lying stellar population probed by the rest-frame optical for 
these galaxies follows the morphology of the Ha emission. 



3.7.1 Radii 

We note that radius is a particularly difficult property to 
measure robustly and consistently with a variety of methods 
used in the literature. We calculate the kinematic plateau 
radius, rp, from disk modeling and the effective radius, 
reft, from surface brightness profile fitting. Th e reg val- 
ues presented are consistent with the sample of iLaw et al.l 
(2009), 0.6 kpc ^ r $; 1.6 kpc, but ar e smaller than most 
othe r samples (e.g. Wright et al.l 20091 : 2.1 kpc ^ r ^ 3.9 
kpc, iForster Schreiber et al.ll2009l : 1.2 kpc < r < 7.5 kpc, 
lEpinat et aT]|2009l : 1.2 kpc ^ r < 5.5 kpc). The haff light 
radius, ri/2, is also estimated directly from the data by de- 
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Figure 9. Example of surface brightness modeling results for 
14s. The black points represent the flux in each spaxel from the 
peak of emission marked with the red plus sign in the Ha image 
of the galaxy in the top right corner. The red circles represent 
binned points, the black solid line is the best-fitting Sersic model, 
and the black dotted line is a fit to the fiux profile of a tip-tilt 
star observed with OSIRIS. 



Table 7. Surface Brightness Fitting 

m n 



refi(kpc) 



Ha Images 



WK0909.02R 


1 n+^-i 


0.61?- 


WK1002.18S 


0.8 ±0.1 


1.9 ±0.4 


WK1002.14S 


0.9laf 


0.9±°- 


WK1003.46S 


1.0 ±0.1 


0.7l°- 


WK1003.52S 


2.0tn 


0.4 ±0.1 


WK0912_16S 


o.5±g;J 


-■-3-0.6 



RCS2 rest-frame UV Images 



WK0912.13R 
WK0909.02R 
WK0909.06R 



TU-r 
1.8 

+0.2 
0.6 



0.6 
1.0: 
0.4 ±0.3 



6.1 ±2.0 



3.4- 



+2.6 



termining the size of the circular aperture needed to enclose 
half the Ha hght from the centre of emission, and is given 
in Table U In some cases ri/2 is the best measure of radius 
as neither retr or rp are able to be measured due to irregular 
morphologies and kinematics. The r-1/2 values are greater on 
average than the r^s , likely due to noise, and are compara- 
ble to the samples of lWright et al.l ([20091) and lEpinat et al.l 
l|2009ll . 



4 DISCUSSION 

The sample presented here was selected from the WiggleZ 
survey for strong UV and nebular luminosities. The lumi- 
nosity, mass, SSFR, and i^shcar are more comparable to the 
z > 2 sub-set of the SINS survey than the samples at the 
same redshift, but more comparable in size to the 2: ~ 2 



sample from iLaw et al.l ((20091). Like other IPS surveys we 
find a range of kinematics and morphologies with high ve- 
locity dispersions on average (((Tmoan) = 98 km s~^), yet the 
uniqueness of this sample is derived from the range of Ha 
morphologies with resolved fcpc-size star-forming regions and 
ultra-compact Ha emission of the single emission galaxies. 

However, IPS samples are all selected differently and 
may not be directly comparable. The SINS survey observes 
galaxies found from a variety of selection tech niques includ- 



ing o ptical color selected 'BX/BM' galaxies (ISteidel et al 



M04|), infrared color selected 'BzK' galaxies ( Daddi et al 



2004), Lyman break selected ga laxies (LBGs ) , and bright 
sub- millimeter galaxies (SMCs). lEpinat et al.l (12009!) select 
their gala xies from the magn itude complete VVDS sur- 
vey, while IWright et al.l ((20091) observe BM galaxies, and 
JLaw et all ( 20091 ) observed BX galaxies. While the WiggleZ 
sample was not selected by any of these methods, the SEDs 
of the sample imply that they could be classified as BM 
galaxies, BzKs and LBGs. Because IPS surveys are observa- 
tionally expensive we are working with low number statistics 
and we cannot follow up mass or luminosity complete sam- 
ples at high redshift. 

The initial classification of this sample was based on 
Ha morphology into three groups; single emission, extended 
emission, and multiple emission. A reclassification is now 
considered based upon the properties discussed in Section 3 
in the context of disk formation at high redshift and a merger 
sequence. However, it is a challenge to distinguish between 
a galaxy merger and the merger of massive sub-clumps em- 
bedded within a rotating disk and similarly to distinguish 
between a merger remnant and a pseudo-bulge formed from 
the collapse of massive sub-clumps. Below we explore models 
of both disk formation and merger events at high redshift. 
We use the kinematic properties, masses, sizes, etc. to dif- 
ferentiate between the dominate formation mechanism for 
the WiggleZ galaxies sample presented here. 

4.1 The Merger Hypothesis 

With separations of 3-10 kpc the Ha complexes may be 
sites of multiple regions of star formation merging from dif- 
ferent parent galaxies. Kinematic quantities need to be uti- 
lized to find support for this hypothesis as deep imaging of 
diffuse continuum emission is unavailable to reveal tidal fea- 
tures characteristic of mergers. We analyze here a scenario 
in which our sample form a merger sequence with pre—, on- 
going, and post— mergers. 

In this scenario the multiple emission galaxies are the 
initial stages of mergers when star formation in separate 
galaxies can still be spatially differentiated. The extended 
emission galaxies are the next stage in the sequence repre- 
senting ongoing merging activity (or pre-mergers with com- 
ponent confusion due to projection effects). The final stage 
in the proposed merger sequence is the single emission galax- 
ies representing merger remnants with high dispersions and 
star formation rates. 

Velocity dispersion maps provide some evidence indi- 
cating that we are seeing separate kinematic components 
merging. The observed peaks in dispersion may correspond 
to centres of individual galaxies or be a result of shock fronts 
between merging systems. Systems 4IS anA 61S show multi- 
ple peaks in velocity dispersion correlated with the brightest 
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regions of emission, indicative that the system is an early- 
stage merger where each galaxy has a peak in velocity dis- 
persion corresponding to an unresolved component. Systems 
13R and 06R show multiple peaks in velocity dispersion cor- 
related with the gas between the brightest regions of emis- 
sion, indicative of the turbulence occurring from the inter- 
action of the gas from the merging galaxies in this scenario 
or the effect of fitting substructure along the line of sight 
with a single Gaussian. 

However, the majority of IFS results classified as merg- 
ers show velocity maps where the velocity does not vary 
smoothly across the source but rather show a quick turnover 
within a resolution element or les s (jNesvadba et al.l 120081 : 
I Wright et al. 2003: Law et al. 2009|). An ideal example, pre- 
sented in iNesvadba et al.l 1 20081 1 . is a merging system ob- 
served as two unresolved "clumps" of [OIII] emission with 
a step in the velocity map between both components. This 
step function, observable with AO resolution, is not present 
in the kinematics presented here. 

The single emission galaxies in this scenario could be 
merger remnants. This comparison is particularly difficult 
as merger remnants produced in simulations show a wide 
range of kinematic and morphological properties with ob- 
servables sensitive to the initial mass ratio of the merger 
IjNaab fc Burkertll2003l ). However, the majority of remnants 
are slow rotators (with v/cr < 1), and have high velocity 
dispersions and minor axis rotat ion (Naab fc Burkcrt 20031 : 
ICox et al.|[200^ : lBois et al.ll2010l) . These properties could be 
used to describe 22S, 16S, US, 18S and 46S (although ma- 
jor and minor axes cannot be distinguished) . These galaxies 
show some indication of rotation but with very low ampli- 
tudes (ushoar ~40 km s~^). But it is difficult to distinguish 
merger remnants and disks. Some simulations even find that 
a low percent age of remnants also have disk-like isophotes 
l|Naab fc Burk crt 2003). 

A fraction of the sample, particularly the lower-mass 
galaxies, are found above the 'main sequence' of star forma- 
tion (Figure [7} , which indicates they may be undergoing a 
burst of star formation typical of a gas-poor or minor merger 
event (JNocskc ct al. 2007}. 

Unfortunately, surface brightness detection limits of 
OSIRIS may introduce a selection effect of detecting only 
the brightest objects in merging systems. This effect can be 
seen in Figure Alb of the Appendix for object 02R where 
an outfiow or companion is seen in the RCS2 image but not 
detected in the Ha OSIRIS image. Classification errors may 
arise as a result of this effect. Deep broadband imaging of 
all sources will help to reveal merging companions. 



4.2 The Clumpy Disk Hypothesis 

We investigate the scenario where the data presented here 
rep resent different s t ages of the clumpy disk models of 
the lElmegreen et al.l (|2008l ) simulations, hereafter EBE08. 
In this context galaxies with multiple regions of emission 
represent the early phases of disk formation when grav- 
itational instabilities result in multiple star-forming com- 
plexes forming under Jeans collapse embedded within a disk 
galaxy. As turbulent speeds decrease relative to the rota- 
tional speed of the disk, the clumps migrate towards the 
centre on a timescale of 1 Gyr, represented here by the 
galaxies with extended emission. Finally, the clumps coa- 



lesce and form a slowly rotating pseudo bulge (|Noguchill999l: 



Dekel et al.l 



Bournaud et al. 200 



ay 

23; 



Elmegreen et al.l 12008. . 

2009bl : iKrumholz fc Dekelll2010 l) , represented here by galax- 
ies with a single spherical emission profile. 

The morphologies of multiple emission systems, 06R, 
13R, and 61S (Figure^ are well matched to the early stages 
of clumpy disk formation (f = 175 — 450 Myr). Systems 06R 
and 13R are well fit by disk models. The clump velocities 
follow the global velocity structure of the galaxy as seen 
in Figure [S] where the major axis of rotation cuts across 
multiple clumps revealing no corresponding deviation in the 
rotation curves jlmmeli et al.)i2004b ). 

The morphologies of 02S, 22S, 4IS and 01 R are either 
consistent with the system being observed at the clump co- 
alescence phase (i = 600 — 800 Myr) or the result of pro- 
jection effects. Systems 02S and 22S, each show two close 
components which can be distinguished when fitting surface 
brightness profiles, while 4IS and OIR have single bright 
extended regions where multiple components cannot be dif- 
ferentiated visually or from surface brightness profiles. 22S, 
41s and OIR are well fit by disk models whereas 02S has 
a fiat velocity map with a peak in dispersion corresponding 
to the brightest peak in Ha emission. The fainter peak of 
emission may be a signature of a galactic outflow resulting 
in confusion in classification. 

Star-forming regions within a disk should have similar 
or consistent metallicities as they are evolving within the 
same conditions as predicted by theor etical models for disks 
at high redshift (|lmmeli et al.ll2004ah . In a major or minor 
merger no similar correlation is expected in the metallicities 
of each galaxy. By integrating the light of each clump, the 
clump metallicities are calculated from Equation [5] Metal- 
licities of clumps in 06R and 61 S are equal within the errors, 
as predicted for the early stages of disk fragmentation. At 
clump coalescence metallicities of individual clumps would 
be consistent as seen in 02S and 22S. Although this sup- 
ports current theories, 12 -I- log(0/H) varies little across the 
sample (A[12 + log(0/H)]~ 0.26). 

Compact single emission galaxies I4S, 16S, 18S, 46S, 
and 52S have morphologies well matched to the later stage 
of simulated formation when the clumps have collapsed to 
form a bulge (i = 1000 Myr). The dispersions of I4S, 46S, 
and 52S are the highest in the sample and are consistent with 
EBE08 bulges, predicted to form by clump coalescence and 
become slow rotators with high turbulence of aioa ~ 130 km 
s~^ in their inner regions. Systems I4S, 16S, 18S all show 
velocity gradients with amplitudes less than 70 km s~^. The 
bulge density profiles in the EBE08 models are well fit by 
Sersic profiles of n ~ 4 in the inner 500 pc and n > 2 at 
greater radii. While we do not resolve the inner 500 pc we 
do resolve to 650 — 750 pc and find n < 2 in all cases for the 
Ha emission. 

Velocity dispersion maps in this sample do not show 
characteristic central peaks in velocity dispersion observed 
in IFS st udies in both local llde Zeeuw et al. 2002^ and high 
redshift (iGenzel et aP I2OO6I : iForster Schreiber et al.l [2OO9I 1 
disk galaxies and therefore cannot be used as an indicator 
for disk classification. As discussed in Section 3.6, this peak 
in velocity dispersion can be attributed to beam-smearing 
across the inner regions of the velocity turnover of the 
galaxy. However, because we are observing only the central 
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region of rotation and have a small PSF as a result of using 
LGSAO, a dispersion peak is not expected. 

Furthermore, if we model the velocity disper- 
sion using the equation for a thick or compact disk 
l|Binnev fc TremaindlJOO^ ): 



/izK(r-) 



(9) 



where h^. is the disk scale height, the expected model 
would give a constant, featureless dispersion map within 
the plateau radius (r ^ rp). This would be consistent with 
galaxies 4^3, I4S, and 13R where a plateau velocity is never 
reached. In this case the observed variations in velocity dis- 
persion are more likely a result of random perturbations 
within the centre of the disk or local perturbations due to 
clumps. 

The simulations of EBE08 and others that model and 
evolve these galaxie s have mainly been developed with 
sticky-particle codes l|Elmegreen et al.|[20oi : lNoguchilll999l ') 
which on ly model the c o ld gas associated with the cold cloud 
medium. Ilmmeli et al.l l|2004al lbl) do include a hot phase to 
their simulation allowing the galaxy models to track stel- 
lar ages and metallicities providing observables that can be 
tested. Furthermore, the initial conditions of such models 
require the disks to have very high velocity dispersions to 
mimic observational results at these redshifts making it dif- 
ficult to distinguish inputs and outputs of the models. More 
realistic models are needed for a more rigorous comparison. 
Regardless of these limitations the data is well-matched to 
the available models. 



5 CONCLUSION 

We favor the clumpy disk model over the merger sequence 
model. We find 7/13 of the sample is well modeled by a 
disk and we find that clumps found in galaxies with multi- 
ple regions of emission do not strongly affect the kinematic 
fi elds. This provides obse rvational evidence for the models 
of lElmegreen et alj l|2008r ) which provide an explanation for 
the irregular morphologies seen at high redshift. 

However, whilst natural to attempt to identify kine- 
matic observations as disks or mergers, as been attempted 
here, this simplified picture of galaxies becomes increasingly 
fallible at high-redshift. More than half of the galaxy kine- 
matics in this sample are well-described by disk models, 
however it is important to be ca utious of definiti ve kinematic 
classifications. For example, Ep inat et al.l (|2010|) provide ev- 
idence that up to 30% of disk galaxies may be misclassified 
as perturbed rotators or complex kinematics with current 
classification schemes. 

It is likely that we are seeing a variety of mecha- 
nisms at work which cannot be broken down into simply 
disks and mergers. Even if disk instability is a dominant 
mechanism in star-formation histories at high redshift, as 
we find here, the disks are still likely to undergo mergers 
in th eir lifetime further complic ating kinematic interpreta- 
tions. |Lenioine:£usserolle_e^_alJ (I201(i), for example, observe 
a gala:xy pair with separation of 12 kpc with one galaxy in 
the pair classified as a disk with velocity gradient of ~91 
km s~^ and the other galaxy classified is compact and likely 
dominated by AGN. It is unclear how to classify this sys- 



tem, which has multiple kinematic components. Further- 
more, With such small number statistics, the varying kine- 
matics across multiple IFS samples could be due to varying 
properties amongst the galaxies selected and instruments 
used, making it difficult to compare results. 

The global properties of the sample, the most luminous 
population, studied at 2 ~ 1.3 with IFS, are in surprising 
agreement with a model of unstable disks. We agree with 
iLaw et al.l i|2009f ) that statistically all high-redshift galaxies 
cannot be mergers. Models show that on average the en- 
hancement of SFR in a merger is only a factor of a few and 
lasts 200 — 400 Myr with peak enhancement lasting for les s 
than 100 Myr (|Tevssier et al.ll2010l : lOi Matteo et al.ll2008l '). 
It has been shown locally that all ultra-luminou s infrared 
galaxies (ULIRGS) are mergers (JDuc et al.l 119971 ) and it is 
therefore possible a high-redshift subset of galaxies are also 
exclusively mergers. However, galaxies on average formed 
their stars more rapidly in the past implying the need for an 
evolving definition of what makes more luminous galaxies 
'extreme'. We therefore stress the importance of the bolo- 
metric luminosity in knowing if the sample presented here 
fall on or above the main sequence of star formation at this 
redshift to aid in making generalizations about galaxy pop- 
ulations at 2 > 1. 

Discerning the dust content will aid in a better under- 
standing of their star-formation histories. Such moderate to 
high mass, gas rich, intensely luminous galaxies are reminis- 
cent of local ULIRGs and z ~ 2 — 3 SMGs. However, critical 
to making this link is to determine their far-infrared and 
sub- mm SEDs, and enabling the measurement of their bolo- 
metric luminosities and dust content. Such quantities will 
inform on whether these galaxies are dust-enshrouded star- 
bursts or dust-free late time proto-galaxies. We will present 
Herschel data in a future paper that will address the dust 
content in this sample. 

We expect to see signatures of outflows in these sys- 
tems as galactic scale winds are expected in galaxies 
forming s tars at > SMp yr ~^and are seen in 75% of 
ULIRGS (|Veilleux et al.ll2005h . Previous studies h ave de- 
tecte d galactic winds with I FS (high-z: t^Sw inban k et al.l 
'200 dl: [Alexander et all I2OI0I: [siiapiro et all l2009l : local: 
ISharp fc Bland-Hawthornl2010l ). The high dispersions to the 
edges of 02S, 52S, and 02R are real and not an inflation 
from lower signal-to-noise data. This increase in dispersion 
towards the edges most likely indicates the presence of galac- 
tic winds which add an underlying broad component of Ha 
(|Veilleux et al.ll2005l ). The influence of winds on this sample 
will be further discussed in a subsequent paper. 

Measuring the properties of the individual clumps 
including size, metallicity, SFR, and velocity dispersion 
to compare to local samples will be important to testing 
models of clumpy disk formation. Only on order ~ 10 
clumps have been resolved with IFS studies using lensing 
and a further ~ 10 resolved or ma rginally resolved within 
the SINS survey JGenzel et al.llioTg '). The properties of the 
10 clumps identified in these systems will be presented in 
our next paper. 

In this paper we: 

• Present kinematic data of 13 star-forming galaxies at 
z ~ 1.3, selected from the WiggleZ Dark Energy Sur- 
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vey, of better or comparable resolution than previous 
studies. 

• Extend the trend of ' Green et al.1 (J2010l ) to higher star 
formation rates and higher dispersions between 1 < z < 
2 (Figure [H). 

• Find the most luminous star-forming galaxies at this 
redshift show a variety of star-forming morphologies 
with a high incidence of star formation in disks. 

• Identify clumps of 1-2 kpc in size at z ~ 1.5 that are 
resolved kinematically with LGSAO IFS without aid of 
lensing. 

• Conclude the clumps in multiple emission galaxies fol- 
low the overall kinematics of the system, and find a high 
incidence of velocity gradients consistent with rotation 
lending support of the clumpy disk scenario of EBE08 
over the major merger scenario. 
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APPENDIX A: 

Al Multiple Emission Galaxies 

A 1.1 WK0912.13R 

WK0912-13R has an extended irregular flux profile with 
faint surface brightness at the centre encircled by a brighter 
ring of emission. This galaxy shows multiple bright regions 
of Ha spatial emission, however due to low signal the in- 
dividual regions are not resolved. The velocity map shows 
an ordered velocity gradient with the "hole" in emission spa- 
tially placed at the central region of the velocity gradient 
and is well fit by a disk model. However, this detection has 
the lowest signal-to-noise of the sample and only a portion 
of the rotation curve is detected. Optical emission is shown 
in Figure [XT] 



A1.2 WK1002.61S 

The morphology of 61S reflects that of "chain galaxies" ob- 
served in the Hubble Ultra Deep Field (UDF) with 2-3 
resolved clumps marked by plus signs in Figure [3] However, 
unlike the clumpy galaxies 06R and 13R in this sample, 6 IS 
has one of the lowest velocity dispersions of this sample of 
o"mean=85 kiu s~^ aud shows little or no velocity separation 
between clumps as shown in Figure IX2l There is a messy ve- 
locity gradient across the length of the "chain" with peaks 
of velocity dispersion spatially aligned with the peaks of the 
eastern most clumps. 



Al. 



WK0909.02R 



WK0909_02R has the most irregular optical morphology of 
the objects with available broadband imaging. The RCS2 
r-band image (0.185" /pixel) shows three distinct regions 
of emission spread over 4 arcseconds, including the central 
brightest region denoted by an x in Figure Ell a lower sur- 
face brightness feature to the NW, and a larger low surface 
brightness feature to the south of the brightest region ex- 
tending ~21 kpc away from the central peak of broadband 
emission. The OSIRIS observations of WK0909.02R has a 
single peak in Hq flux with a faint tail of emission to the NW 
coincident with the peak of r-band emission and faint NW 
feature respectively. [Nil] A A 6584 is detected with peak flux 
coincident with peak Ha flux. Long-slit spectroscopy with 
LRIS reveals the southern low surface brightness component 
is at the same redshift as the brighter central emission. 



A 1.4 WK0909.06R 

WK0909_06R has three resolved knots of Ha emission. In 
velocity space these peaks correspond to a velocity range 
of approximately 320 km s^^ shown in Figure IA2I This 
large change in velocity over ~10 kpc leads to broad veloc- 
ity width of the Ha emission line between flux peaks where 
a signature of two Gaussians is seen. Peak [Nil] emission 
corresponds with the brightest peak of Ha and is not de- 
tected elsewhere in the galaxy. Faint [SII] emission is de- 
tected in the integrated spectrum but does not have high 
enough signal-to-noise to be mapped spatially. The RCS2 
image in Figure ETI reveals a possible optical companion 1.8" 
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Figure Al. RCS2 r-band and OSIRIS Ha images for 02R, OIR, 
06R, and 13R, the only objects with resolved broadband imaging. 
The green dot-dash lines on the RCS2 image of 02R mark the 
position of the LRIS long-slit. LRIS spectroscopy confirms the 
low surface brightness region below the central emission is at the 
same redshift as the central emission detected by OSIRIS. 
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to the south-west (15 kpc). Future long-slit spectroscopy will 
determine if the galaxies are under-going a merger. 

A2 Extended Emission Galaxies 

A2.1 WK0905.22S 

The flux proflle of WK0905.22S is also best fit by two com- 
ponents with the the primary peak of velocity dispersion 
occurring at the boundary of the two systems indicative of 
a line-of-sight merger, however the velocity map shows evi- 
dence of rotation and is well fit by a disk. 



A3. 2 WK1002.US 

WK1002-14S has a single compact region of resolved Ha 
emission showing a linear velocity shear, with some varia- 
tions in the velocity field at the edges of the galaxy detec- 
tion. As with ^^5 there exists a degeneracy between Vp and 
Vp due to the linear nature of the velocity profile. The fiat 
dispersion map can be explained by the lack of turn-over in 
the velocity profile. This galaxy is one of the highest dis- 

-^130 km s~^ and stellar mass 



persion galaxies with a^ 



of M.=3.7x 10'" M0 
with Ha emission. 



[NII]AA6583.5 is detected coincident 



A2.2 WK0912.01R 

WK0912_01R has one extended peak of Ha flux with a 
fainter surrounding halo. The velocity map shows a strong 
velocity gradient of Ushoar ~ 200 km s~^ and a near linear ro- 
tation curve. It is well-fit by a disk with optical morphology 
consistent with being an edge-on system as seen by the el- 
lipticity of the rest-frame UV image in Figure IXT] [NH] was 
not detected in the OSIRIS observation of this galaxy, which 
is partly due to its short integration time, an upper limit is 
given in Table [S] 



A2.3 WK1002MS 

The Ha morphology oi ^IS reveals two peaks of emission 
connected by a dip in the Ha surface brightness at the cen- 
tral region of the detection. The rotation profile of ^iS is 
linear (t;shcar=114 km s~^) with a featureless velocity dis- 
persion; as a result there exist a strong degeneracy between 
Tp and Vp when fitting a disk profile. This indicates that 
41s is well-fit by a disk but as we only detect the inner part 
of rotation, the kinematic extent cannot be constrained. It 
follows that the dispersion map is fiat and featureless as the 
turnover has not been reached. 



A3. 3 WK1002U6S 

WKIOO2.46S is the most compact ( ri/2(Ha)~1.8 kpc ) re- 
solved galaxy. The stellar mass oi 46S is 7 x 10^ M© and the 
galaxy has a relatively low amcan— 85 km s~^. 



A3.4 WK1002.52S 

52S is only marginally resolved and an outlier in many re- 
spects including: highest velocity dispersion (o"inean=125 km 
s~^), highest stellar mass (M*=4.6 x lO'^^ M©), most com- 
pact (ri/2 ~1.8 kpc). Given the high stellar mass it is more 
likely that 52S is a merger remnant rather than a newly 
forming bulge. 



A3. 5 WK0912A6S 

WK0912-16S is observed to have a single resolved peak in 
Ha flux. [NH] is detected with peak flux coincident with the 
peak Ha flux. 16S has one of the lowest velocity dispersions 
of the sample at (Tmcan~ 90 km s~', and does not show a 
coherent velocity gradient. It also has the lowest stellar mass 
of the sample of 6.3 x 10^ M© and the flattest flux profile 
with nHa=0.5J 



-+0.3 

-0.1- 



A2.4 WK0809.02S 

The 2-D fiux profile of WK0809.02S is best fit by two emis- 
sion components. The peak in velocity dispersion corre- 
sponds with the brightest peak in flux, however no clear 
velocity gradient is seen and the object has the lowest ve- 
locity shear in the sample (wahear =36 km s~^). With J band 
photometry {K not available) the stellar mass is estimated 
at log(M,) = 11.1 ± 0.3 Mq, one of the most massive of the 
sample. A tentative IRAM detection of the C0(2-l) emission 
line exists at 500 km s~' from the Zgys of Ha. The CO(2-l) 
line is formally 2.5o- per channel, but spans >3 channels, 
so is ~5(7 integrated over the line, possibly indicative of an 
outfiow/infiow signature. This result is outside of the scope 
of this analysis and will be addressed in a subsequent paper. 



A3 Single Emission Galaxies 

A3.1 WK1002A8S 

WK1002-18S is a compact resolved single emission galaxy 
with nHa=0.8 and shows a near-linear velocity gradient. The 
stellar mass of 18S is 8 x lO^M0and has a velocity dispersion 

0-mcan= 102 km S~^. 
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Figure A2. Integrated spectra from individual clumps for all 
multiple emission systems. Clumps and their corresponding spec- 
tra are labeled in Ha images found in the upper right corner. 
The spectra are the integrated spectra of emission within 2 r^g 
of star-forming regions that are resolved or marginally resolved. 



